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ABSTRACT 
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average high school student who is enrolled in a Pascal course or who 
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Drawlinecieneral, DrawCircleGeneral, plotcuces, ConvertToPolar, and 
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Introduction 

A Potpourri of Pascal Programs is a collection of Pascal programs that were developed for an 
NSF high school teacher's workshop which was held in the summer of 1986. Some of the programs 
are standard examples in computer science and no source is indicated for these programs. Other pro- 
grams were taken from published literature and, in these cases, the original source is indicated. Most of 
the published programs were written in BASIC and translated by us into Pascal. The remaining pro- 
grams represent original work by either the workshop staff or by the woikshop participants. The 
appropriate author is indicated for these programs. Some of the original programs utilize graphics pro- 
cedures which are given in the appendix and were written by Dr. Ziegler. 

The programs can be used as a means of extending or enriching textbook material in either high 
school mathematics or Pascal courses. Some suggested uses are: 

(1) teacher demonstrations in mathematics classes 

(2) programs for student use in mathematics classes 

(3) student assignments in a mathematics or Pascal course. 

The topics which are addressed in the Pascal programs come from a variety of mathematical 
areas: algebra, trigonometry, discrete mathematics, number theory, mathematics^ modeling and numeri- 
cal algorithms. Generally speaking, the programs are quite short and the prograi.iming level is appropri- 
ate for the average high school student who is enrolled in a Pascal course or who has completed such a 
course. 



ALGEBRA PROGRAMS 

These programs deal with topics in high school algebra courses. 



ALGEBRA PROGRAM 1: Binomial Squares. 
AUTHOR: Lynn R. Ziegler 

OBJECT To dispel an incorrect notion which students often hold about a binomial 
square. 

e.g. (x+2)2 ;t + 4 but (x+2)2 = + 4x 44 



Listing of Binomial Squares Program 

program binomialSquare( input, output ); 

var x: integer; 

begin {binomialSquare} 
wiite!n(lst,' X (X+2r2 X^2-h4 r2-H4X-H4'); 

writeln(lst,' ~ 

for x:« -5 to 5 do 
writein(lst,x:2,sqr(x+2):9, 

(sqr(xM):7, 
(sqr(xM*x-H4):10); 
end. {binomialSquare} 
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Sample Output 



X 


(X+2)*2 X*2-h4 


X*2+4X+4 


-- 

•5 


9 


29 




9 


•4 


4 


20 


4 


•3 


1 


13 


1 


-2 


0 


8 


0 


•1 


1 


5 


1 


0 


4 


4 


4 


1 


9 


5 


9 


2 


16 


8 


16 


3 


25 


13 


25 


4 


36 


20 


36 


5 


49 


29 


49 



ALGEBRA PROGRAM 2: Mixture 



AUTHOR: Randy Odendahl (based on an idea in the 1984 NCTM Yearbook, page 
187) 



Object To compute the cost of a mixture resulting from various proportions of two 
different coffees. (This is an extension of the prob! ^p. usually solved in algebra). 



Typical problem: If Coffee Type I costs $4.00/lb and Coffee Type II costs $5.40/lb, 
how many ounces of Coffee I and how many ounces of Coffee II should you combine 
to make a mixture selling for $4.35/lb? 



Mixture] 
. 16oz J 



16 



16 



16 



4.00 + 86.40 - 5.40 = 69.60 



-1.40 jc = -16.80 



-16.80 

X = — - = 12 02 

-1.40 
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Natural Extension: What would be the cost/lb for various propcxtions of the two 
different coffees? 



4.00 ;c + 5.40(I6--;c) = 16 MIX 



4,00 X 5.40(I6-;c) 



Listing of Mixture Program 

program mixture; 

{mixture computes cost of mixture resulting from 
various proportions of two different coffees.} 

var coffeel: integer; 

begin 

wriieln(lst,'COFFEE I COFFEE II Cost/pound'); 

wriieln(lst,'- — '); 

for coffee 1> 0 to 16 do 

writf!in(lst,coffee 1 :5,( l6-coffee I): 1 0, 
((4.0/I6.0)*coffeel+ 
(5.4/l6.0Hl6-coffeel)):l3:2); 

end. 
Sample Output 



COFFEE I 


COFFEE II 


Cost/pound 


0 


16 


5.40 


1 


15 


5.31 


2 


14 


5.22 


3 


13 


5.14 


4 


12 


5.05 


5 


11 


4.96 


6 


10 


4.87 


7 


9 


4.79 


8 


8 


4.70 


9 


7 


4.61 


10 


6 


4.52 


11 


5 


4.44 


12 


4 


4.35 


13 


3 


4.26 


14 


2 


4.17 


15 


1 


4.09 


16 


0 


4.00 
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ALGEBRA PROGRAM 3: Rectangle 

AUTHOR: Bruce Carlson - Dciiar Bay High School, Dollar Bay, Michigan 49922 
OBJECT Comi>ute the width and an^a of a rectangle given the perimeter and length. 
Listing of Rectangle Program 

program rectangle(input,output); 

{This program is designed to tell a person what the width and area of a 
rectangle would be given the perimeter and length.} 

var length,width,perimeter,area:real; 

procedure findwidth(perimeter,length:real); 
begin {findwidth} 

width>(perimeter/2)-length; 
end; {findwidth} 

procedure findaiea (length,width:real); 
begin {findarea} 

area>length*width; 
end; {findarea} 

begin {rectangle}; 
clrscr; 

writeln ('This exercise is designed to give you the width and the area of 
'a rectangle once you have decided upon a perimeter and length.'); 

writeln; 

writelnCTry to get the largest area for a given perimeter by changing the length.'); 
writeln; 

writeln('What is the perimeter that you would like?) It must be positive.'); 
readln (perimeter); 

writeln('What is the length that you would like?'); 

readin(length); 

writeln; 

if ((length<(perimeter/2.0)) and (length>0.0)) then begin 
findwidth(perimeter, length); 
writelnCThe width of your rectangle is ',width:20:4); 
writeln; 

findarea(length,width); 

writeln('The area of your rectangle is length times width or', area:20:4); 
writeln; 

writeInC Write the perimeter, length, width and area on a sheet of paper.'); 
writeln; 

writelnCTry different lengths without changing the perimeter.'); 
writeln; 

writeln('What happens to your area as your length and width become', 
' closer to each other?'); 

writeln; 

writelnCPress the run key (R) before trying a new length.'); 
end {if} 
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else begin 

writelnCTHINK! Your head is not just a hair fann«'); 
writeto('Did you choose a POSITIVE PERIMETCR?'); 
wriiebCAlso, you must choose a positive leagth that is less than\ 

' one-half of the 
wriielnCperimeter that you choose. Do you know why?'); 
wriieln; 

writelni'Now press the run key(R) and try again.'); 
end;{if-then-else} 
end.{rectangle} 



Sample Output 

This exercise is designed to give you the width and the area of a rectangle once 
you have decided upon a perimeter and length. 

Try to get the largest area for a given perimeter by changing thelength. 

What is the perimeter that you would like?) It must be positive. 
14 

What is the length that you would like? 
8 

THINK! Your head is not just a hair farm. 

Did you choose a POSITIVE PERIMETER? 

Also, you must choose a positive length that is less than one-half of the 

perimeter that you choose. Do you know why? 

Now press the run key(R) and try again. 

>r 

This exercise is designed to give you the width and the area of a rectangle once 
you have decided upon a perimeter and length. 

Try to get the largest area for a given perimeter by changing thelength. 

What is the perimeter that you would like?) It must be positive. 
14 

What is the length that you wou!d like? 
4 

The width of your rectangle is 3.3000 

The area of your rectangle is length times width or 12.0000 

Write the perimeter, length, width and area on a sheet of paper. 

Try different lengths without changing the perimeter. 

What happens to your area as your length and width become closer to each other? 



Press the run key (R) before trying a new length. 
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ALGEBRA PROGRAM 4: Shipping 

AUTHOR: J. Bugni, L'Anse High School, L'Anse, Michigan 49946 

OBJECT Determines if UPS will ship your package given its width, height, and 
depth. 

Listing of Shipping Program 

{ Shipping program by J. Bugni, L'Anse High School } 
program shipping (input,output); 

{This propogram determines if UPS will ship your package} 
var w,h,d,total:real; 

{function uprmeasure finds the sum of the girth and length of a package.} 
function upsmeasure (w,h,d:real): reM; 

{function upsmeasure will find the guth of a package and add to this 
number the length of the package} 

var I,gl,g2,g3,girth:rcal; 

begin {upsmeasure} 
l:.w; 

if h>l then l>h; 
if d>l then l:»d; 

gl:-2*w+2*h; 
g2:-2*h+2*d; 
g3>2*w+2*d; 

girth>gl; 

if g2<girth then girth>g2; 
if g3<girth then girth:«g3; 

upsmeasure:«girth+l; 
end; (upsmeasure/ 

begin {shipping} 

write ('enter the width of your package > '); 

readln (w); 

writeln; 

writeln; 

write ('enter the height of your package > '); 

readln (h); 

writeln; 

writeln; 

write('enter the depth of your package > '); 
readln (d); 
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writein; 
writein; 

total:«upsmeasure (w,h,d); 

if total>108 then write ('sorry your package is unaccepuble') 
else begin 

write (*we are happy to send your package*); 

writein; 

writein; 

writein ('hinl:stand your package up so that the longest side is'); 
writeln('the vertical side, and put the label on the top!!!'); 
end; {else} 
end. {shippiiig} 



Sample Ouiput 



Running 

enter the width of your package > S 



enter the height of your package > 12 



enter the depth of your package > 9 



we are happy to send your package 

hint:stand your package up so that the longest side is 
the vertical side, and put the label on the top!!! 

Running 

enter the width of your package > 43 



enter the height of your package > 43 



enter the depth of your package > 43 



sorry youi* package is unaccentable 
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TRIGONOMETRY PROGRAMS 



TRIGONOMETRY PROGRAM 1: Variable Sine Curve 

AUTHOR: W. Galfhey, Hancock Public Schocb, Hancock, Michigan 49930 

OBJECT To change amplitude and period and graphically display a phase shifted sine 
curve. 

Listing of Variable Sine Curve Program 
program variablesinecurve; 

{NSF CLASS PROJECT— -W. GAFFNEY} 

{This program allows you to change the ampliaide and period 

and cause a phase shift of the standard sine curve} 

var ij:integer; 
k,r,a,b,c:real; 

{Sigraphics} 
{$igraphics.two} 

begin 

^^ritelnCTHIS PROGRAM ALLOWS YOU TO CHANGE THE AMPLITUDE'. 

' AND PERIOD'); 
writelnC ALONG WITH A PHASE SHJFT OF THE STANDARD'. 

' SINE CURVE'); 
writeln('Y-aSIN(bX+c)'); 

writelnCTHE VALUE OF a WELL DETERMINE THE AI^IPLITUDE'); 
writelnCTHE VALUE OF b WILL DETERMINE THE PERIOD'); 
writelnCTHE VALUE OF c WILL CAUSE A PHASE SHIFT'); 
writelnCTHE STANDARD CURVE HAS a-1. b-1. c-0 AND IS'. 
' SHOWN IN RED.'); 

writelnCINPUT THE VALUE OF a'); 
readln(a); 

if (a>2.5) or (a<-2.5) then begin 

wrileln( 'value of a must be in the range of -2.5 to 2.5'); 

readln(a); 
end;fif} 

writelnCINPUT THE VALUE OF b'); 
readln(b); 

writelnCINPUT THE VALUE OF c'); 
readln(c); 
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graphicson; 
cleargraphic<; 

{plot tha axes} 

for i>0 to xmax do plotgeneral (1,100,0,319,0,199,1); 
for i>0 to ymax do plotgeneral (160,i,0,3 19,0, 100,1); 

for i>0 to 639 do begin 
n-i*pi/180; 

{plot the standard curve} 

k:-^0*sin(1.12*r)+100; 
plotgeneral(iJcO,639,0, 199,2); 

{plot the variable curve} 

k>a*40*sin((b* 1.12*r)-H:)+100; 
plotgeneral (i,k,0,639,0, 199,3); 
end; 

delay(lOOOO); 
graphicsoff; 
end. 



Sample Output 
Running 

THIS PROGRAM ALLOWS YOU TO CHANGE THE AMPLITUDE AND PERIOD 
ALONG WITH A PHASE SHIFT OF THE STANDARD SDJE CURVE 
Y-aSIN(bX-K:) 

THE VALUE OF a Wtti DETERMINE THE AMPLITUDE 

TOE VALUE OF b Wtti DETERMINE THE PERIOD 

THE VALUE OF c Wtti CAUSE A PHASE SHIFT 

TOE STANDARD CURVE HAS a-1, b-1, c-O AND IS SHOWN IN RED. 

INPUT THE VALUE OF a 

2 

INPUT THE VALUE OF b 
1 

INPUT THE VALUE OF c 
5 
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DISCRETE MATHEMATICS 



DISCRETE MATHEMATICS PROGRAM 1: Factorial and Recunive Factorial. 



AUTHOR: Randy CJendahl 



OBJECT To contrast the use of an iterative algorithm with the use of a recursive 
algoriu\m for generating n!. 



Listing of nonrfcursive factorial: 

program factorial(input,output); 
var t: integer; 

n: real; 
begin {factorial} 

writcln(lst,*What number would you like the factorial of?*); 

readln(n); 

write(lst,n)und(n):l/ factorial is '); 
for t:« n>und(n) downto 2 Jo 

n:- n*(i.l); 
writeln(lst,ix20:0); 
end. {factorial} 



Listing of recursive factorial: 

program recursiveFactorial(input,output); 
var 
n:integer; 

function factorial( n: integer ): real; 
begin {factorial} 

if n«0 then facu»ial:« 1 

else factonal:- n*factorial(n-i) 
end; {factorial} 

begin {execute recursiveFactorial} 
writeln(Ist,'What number would you like *. 
'the factorial of?*); 

readln(n); 

writeln(output,n:l,* factorial is ',factorial(n):20:0); 
end. . cecute recursiveFactorial} 



twisting of sample output: 

What number would you like the factorial of? 
17 

17 factorial is 568742809600 
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DISCRETE MATHEMATICS PROGRAM 2: Fibonacci & Recursive Fibonacci 
AUTHOR: Randy Odendahl 

OBJECT To contrast the use of an iterative algorithm with the use of a recursive 
algorithm for generating Fibonacci numbers. 

Listing of nonrecursive Fibonacci: 

program Fibonacci; 

{Fibonacci prints out the first 20 Fibonacci numbers.} 

var al,a2,a3,i: integer; 
begin {Fibonacci} 
al:- 1; a2:- 1; 
writeln(kt,' i ith Fibonacci'); 

writeln(Ist,' - '); 

writeb(lst,l:2,al:10); 
writeln(Ist^:2,a2:10); 
for i:- 3 to 20 do begin 
a3:- al +a2; 
writeln(lst,i:2,a3:10); 
al:- a2; 
a2:- a3; 
end {for}; 
end {Fibonacci}. 

Listing of recursive Fibonacci: 

program recursiveFibonacci; 
var i: integer; 

function Fibonacci( n: integer ): integer; 
begin {Fibonacci} 

if (n«l) or (n«2) then Fibonacci:* 1 

else Fibonacci:- Fibonacci(n-l)+Fibonacci(n-2); 
e-'d; {Fibonacci} 

begin {recursiveFibonacci} 

writeln(lst,' i ith Fibonacci'); 

writeln(lst,' - '); 

for i:- 1 to 20 do writein(lst,i:2J^ibonacci(i):10); 
end {recursiveFibonacci} 

Listing of sample output: 

i ith Fibonacci 



I 
2 
3 
4 
5 



I 
I 
2 
3 
5 



15 
16 
17 
18 
19 
20 



610 
987 
1597 
2584 
4181 
6765 
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DISCRETE MATHEMATICS PROGRAM 3: Computing powers. 



AUTHOR: Lynn R. Ziegler 



OBJECT To show a different type of recunive function. (In this case, one useful for 
conq)uting int^er powers of real numbers • x". It works by using the observation 
that x2* = (x*/andx^*> = x (x*)l) 



Listing of recursivePower: 

program recursivePdwer(input,output); 

var x:re3!; 
n:integer; 

function power(x:real; n:integer) : real; 
var temp:real; 
begin {power} 

if n-0 then power- 1 
else begin 

ten^:-power(x,n div 2); 
if (n mod 2) - 0 then power«sqr(temp) 
else power>sqr(temp)*x 
end {if then else} 
end; {power} 

begin {recursivePower} 
writeln('Enter x and n '); 
readln(x,n); 
Wiitttin(x:8:2,' to the', 
n:3,'th power is ', 
power(x,n):18:I) 
end {recursivePower} 



Sample Output: 

Enter x and n 
2.0 50 

2.00 to the SOth power is II2S899906842620.0 



ERLC 



15 



A Potpourri of Pascal Programs 



PageU 



NUMBER THEORY 

NUMBER THEORY PROGRAM 1: Printmod 
AUTHOR: Lynn R. Ziegler 



OBJECT To print the fint 100 counting numbers modulo the integer input by the 
user. 



Listing of Printmod: 

program printmod(input,output); 

{printmods ou^ts the first 100 counting numbers modulo die 
integer input by the program's user.} 
var ij: integer; 
begin 

writeln(lst,'The first 100 coundng numbers \ 

'in modular arithmetic'); 
writeln(lst,' 

writeln(lst,'Please enter an integer for the modulus:'); 

readln(j); 

writeln(lst); 

writelnOst,' n ',' n mod 'j:l); 

wnteln(lst,' - ',' '); 

fc i> 1 to 100 do 
writeln(Ist,i:4,(i mod j):4) 

end. 



^ar'iip': Output: 

Ik . ^ rs; .X) counting numbers in modular arithmetic 
V'•-^.^c enter an integer for the modulus: 

n n mod 6 

1 1 93 3 

2 2 94 4 

3 3 95 5 
-1 4 96 0 

5 5 97 1 

6 0 98 2 

7 1 99 3 

8 2 100 4 
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NUMBER THEORY PROGRAM 2: Wondrous 

SOURCE: [MiUigan], p. 99 translated to Pascal by Randy Odendahl 

OBJECT To determine if a given number is "wondrous". Begin with a whole 
number. If odd, multiply by 3 and add 1. If even divide by 2. If repeated applica- 
don of this proctdurt yields 1 then the original number is called "wondrous". If the 
given number is not wondrous then the program never halts. 

For example, if we stait with 102: 

102/2-61 61*3+1-184 184/2-92 92/2-46 46^2-23 

23*3+1-70 

4/2-2 2/2-1 (after 20 steps) 

Listing of Wondrous 

program wondrou$(input,output); 
var iterationCount»candidate,temp: integer; 
done: bodean; 

begin 
done:- false; 
while not done do begin 

writehi(Ist,*PIease enter number to be tested,*, 
* or enter 0 to stop.'); 

read(candidate); 
if candidateoO then begin 
iteradonCount:- 0; 
temp:- candidate; 
while tempo 1 do begin 
if odd(ten4>) then begin 
iteradonCounc- iterationCount +1; 
temp:- tttnp*3 +1; 

end {if odd}; 

iteradonCount:- iteradonCount +1; 
temp:- temp div 2; 

end {while}; 

writebi(Ist,iterationCount,* iterations were ', 

'needed to make ',candidate,' wondrous.'); 
end {if candidate} 
else done:- true; 
end {while not done}; 
end 

Listing of sample output: 

Please enter number to be tested or enter 0 to stop. 
102 

20 iterations were needed to make 102 v^ondrous. 
Please enter number to be tested or enter 0 to stop. 
0 
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MATHEMATICAL MODELLING 

MATHEMATICAL MODELLING PROGRAM 1: ExpGrowth 



AUTHORS: Beverly Ginunestad and Randy Odendahl 



OBJECT Model growth of the US population from 1790-1970 using the exponential 
function P{t) = 3.929 < Population estimates are printed at ten-year intervals. 



Listing of ExpGrowth: 

I»iogram expGrowth(input,output); 

const baseYear> 1790; 

stopYear- 1970; 
var elapsedTime : integer; 
begin 

vmteln(lst,'YEAR POPULATION(mmions)'); 

elapsedTime:> 0; 

repeat 

writeln(lst,baseYear+el^sedTinie:4, 

3.929*exp(0.029655*elapsedTime):15:3); 
elapsedTime:> elapsedTime +10; 
until ( (baseYear -felapsedlime) > stopYear ); 
end. 



Sample Output: 

YEAR POPULATION(nullions) 

1790 3.929 

1800 S.28S 

1810 7.110 

1820 9.S64 

1830 12.866 

1840 17.307 

18S0 23.282 

1860 31.319 

1870 42.131 

1880 56.675 

1890 76.240 

1900 102.559 

1910 137.963 

1920 185.589 

^30 249.656 

1940 335.840 

1950 451.775 

1960 607.733 

1970 817.528 
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MATHEMATICAL MODELLING PROGRAM 2: Prey-Predator 



SOURCE: [Wapner], pp. 137-8 translated into Turbo Pascal by Randy Odendahl. 



OBJECT Present a mathematical model describing the population of rabbits (r) and 
the population of wolves (w) using difference equations. (See The Mathematics 
Teacher, February 19H 137-138.) 

Note: Foi the differential equaions approach see Kemeny and Snell's Mathematical 
Models in the Social Sciences. Blaisdell Publishing Co» 1962. 

In the absence of wolves the population rabbbits would grow at a rate propcxtional 
to its size (At = a r, a > 0). In the absence of rabbits* the population of wolves 
would die at a rate proportional to its size (Aw ^-d w, d >0), When the two popu- 
lations interact, the population of rabbits will be decreased by a term pnqxKtional to 
the number of kills and the population of wolves will be increased by a term propw- 
tional to the number of kills, llie number of kills will vary jointly as r and w. 
Ar=flr - b r w 

Aw=c rw - dw "^^'^ <^*b.c.d > 0 



Listing of preypredator: 

program preypredator(input,output); 

{delta r • ar-brw, delta w« crw-dw. 
the values of a»b,c,and d given in the dau 
statement below will give an equilibrium 
point at 300 rabbits and 200 wolves} 

const a»0.04;b»0.0002;c«0.0001;d«0.03; 
var nw:real; 

t: integer; 
begin {pitypredator} 
writebi(lst,'Enter initial number of rabbits', 

' and wolves'); 
read(inputtr,w); 

writeln(lst,' 25 year period rabbits wolves'); 
writeln(lst); 

writeb(Ist,round(r):7, round(w):7); 
for t> 1 to 1000 do begin 

r:« r +(a*r -b*r*w); 

w:« w + (c*r*w -d*w); 

if t mod 25 « 0 then 
writeln(lst,(t div 25):8,round(r):12,round(w):7); 

end; 

end. {preypredator} 
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Sample Output: 

Enter initial number of rabbits and wolves 
350 120 

25 year period rabbits wolves 



1 


350 


120 


22 


211 


130 


*> 

40 


482 


163 


23 


313 


117 


3 


464 


262 


24 


451 


145 


4 


284 


314 


25 


494 


232 


S 


179 


256 


26 


333 


313 


6 


162 


183 


27 


197 


* 278 


7 


201 


134 


28 


161 


201 


8 


294 


117 


29 


185 


145 


9 


431 


137 


30 


263 


119 


10 


502 


216 


31 


392 


127 


11 


362 


308 


32 


502 


188 


12 


210 


288 


33 


415 


290 


13 


162 


212 


34 


241 


304 


14 


178 


151 


35 


167 


233 


15 


248 


120 


36 


168 


165 


16 


371 


123 


37 


223 


126 


17 


494 


174 


38 


332 


118 


18 


441 


277 


39 


469 


154 


19 


261 


310 


40 


481 


248 


20 


172 


245 


41 


306 


314 


21 


165 


173 









MATHEMATICAL MODEUNG PROGRAM 3: Buffalo Simulation 



SOURCE: Dwayne ChanneeU and Christian Hirsch, "Computer Methods for Problem 
Solving in Secondary School Mathematics", 1984 SCTM Yearbook, pp. 178-181. (See 
Bibliography for a full reference for the yearbook.) 

Translated to Pascal by Randy Odendahl. 



OBJECT The object of this program is to develop a mathematical equation which can 
predict buffalo herd size over the next ten years, based upon knowledge of the current 
number of adult males, adult females, male calves, and female calves. 

Note: It was desired to have a harvesting policy which would not actually endanger 
the buffalo population and this program was developed to assist in the formulation of 
such a policy. 
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Listing of Buffalo Simulation 

program BufraioSimulation(input,output); 
var 

{SIMULATION VARIABLES} 
adulcMales,adultFemaies^eCaives,femaleCalves. 
babyMale$,babyFemales,yearlingMales»yearlingFcmaies, 
herdSize: real; 

yean integer; 

{TEMPORARY STORAGE DURING CALCULATIONS} 
tAdultMaies,tAdultFemales,tBabyMaies,tBabyFemales, 
tYearlingMaies^tYearlingFemaies: real; 

{FORMATTING VARIABLES} 
i: integer; 

begin 

writeln(lst,'Enter number of adult males'); 
readln(adultMales); 
writeln(lsuound(adultMales):I); 
writeln(lst,'Enter number of adult females'); 
readln(adultFemales); 
writeln(lsuound(adultFemales): I); 
writeln(lst,'Enter number of male calves'); 
readln(maleCalves^ 
writelnOsMoundCinaleCalves): I); 
writeln(lst,'Enter number of female calves'); 
readln(femaleCalves); 
writeln(lsUound(femaleCalves): I); 

{MODEL ASSUMES TWO-THIRDS OF CALVES ARE NEWBORN 
ONE-THIRD I YEAR OLD} 
babyMales:- 2.0^3.0*maleCalves; 
yearlingMales:- maleQlves-babyMales; 
babyFemales:- 2.0/3.0*femaleCalves; 
yearlingFemales> femaleCalves-babyFemales; 

writelnOst); writelnflst.' BUFFALO POPULATION DISTRIBUTION'); 
writeln(lst.'HERD':6,' ADULT': 10,' ADULT': 10,'MALE': 10,'FEMALE': 10); 
writelnast,'SIZE':6,'MALES':10,'FEMALES':lO,'CALVES':lO/CALVES':IO); 
for i:« I to 50 do 

write(lsu'-'); 
writeln(lst); 
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{BEGIN SIMULATION} 
for yean* 1 to U do begin 
{COMPUTE HERD DISHUBUTTON FOR EACH YEAR} 
herdSize:* adultMales+adultFemales+babyMales+babyFeinales+ 

yearlingMales+yearlingFemales; 
writeln(Isuound(herdSize' :6, 

round(adultMales):10, round(aduitFemales):10, 
round(babyMales+yearlingMales): 10, 
n)und(babyF^inaIes+yearlingFeinale$):10); 

{CALCULATE NEW DENSITIES} 
tAdultMales:- adultMales; 
tAdultPemales:* adultPemales; 
tBabyMales> babyMales; 
tBabyFemales:* babyFemales; 
tYearlingMales:- yeariingMales; 
tYearlingFenuUes> yearlingPemales; 
adultMales:- 0.9*tAdultMales -i^.6*tYearlingMales -1000; 
adultFemales> 0.9*tAdulcFdnales 40.6*tYearlingFeinales; 
babyMales:- 0.48*tAdultFeinales; 
babyFemales:- 0.42*tAdiiltFeinales; 
yeariingMales:- 0.5*tBabyMales; 
yearlingPemales:- 0.5*tBabyFemales; 
end {for}; 
end {program buffaloSimulation} 



Sample output: 

Enter number of 9/iult males 
10400 

Enter number of adult females 
9100 

Enter number of male calves 
3380 

Enter number of female calves 
3120 



BUFFALO POPULATION DISTRIBUTION 
HERD ADULT ADULT MALE FEMALE 
SIZE MALES FEMALES CALVES CALVES 



26000 


10400 


9100 


3380 


3120 


28207 


9036 


8814 


5495 


4862 


28393 


7808 


8557 


6415 


5613 


27853 


7338 


8848 


6223 


5445 


27760 


6873 


9073 


6300 


5513 


27810 


6418 


9244 


6479 


5669 


27888 


6050 


9435 


6615 


5788 


28037 


5752 


9634 


6747 


5904 


28260 


5508 


9836 


6889 


6028 


28544 


5316 


10041 


7033 


6154 


28885 


5171 


10251 


7180 


6283 
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MATHEMACIAL MODELUNG PROGRAM 4: Falling Bodies 



AUTHOR: Ernest Mattson, Ironwood Area Schools, Ironwood. Michigan 49938 



OBJECT Model the behavior of falling bodies according to Newton's law. 



Listing of Falling Bodies 



PROGRAM FALLINGBODIES (INPUT.OUTPUT); 

{THIS PROGRAM WILL FIND THE DISTANCE THAT A ^TIEELY 
FALLING BODY WILL FALL IN THE TIME PERIOD YOU 
SPECIFY. IT WILL MAKE A CHART OF THE COMPUTED 
VALUES AND DRAW THEIR GRAPH. YOU MUST ENTER THE 
BEGINNING AND ENDING TIMES AS REAL NUMBERS EQUAL 
TO OR GREATER THAN ZERO. WHEN YCU ENTER A 
BEGINNING NUMBER GREATER THAN ZERO THE COMPUTED 
VALUES AND THE GRAPH WILL SHOW THE INTERVAL OF TIME 
THAT YOU ENTERED AND THE CORRESPONDING DISTATICES.} 

CONST G-9.8; 

VAR S:REAL; TREAL; TiREAL; 
LAST:REAL; 
DELTATREAL; 

{$IGRAPHICS.PAS} 
{$IGRAPHICS.TWO} 

PROCEDURE PRCDRW (XO,YO,X1,Y1,COLOR:INTEGER): 

VAR SLOPE:REAL; 
X,Y:INTEGER; 

BEGIN 
IF (XloXO) THEN BEGIN 
SLOPE>(Y1-Y0)/(X1-X0); 
IF (X1>X0)THEN BEGIN 
FOR X>X1 DOWNTO XO DO BEGIN 

PLOTPOINT (XJ10UND(SL0PE*(X-X1 )+Yl ),COLOR); 
END; (♦FOR*) 
END (♦IF*) 
ELSE IF (X0>X1)THEN BEGIN 
FOR X>XO DOWNTO XI DO BEGIN 

PLOTPOINT (XJIOUND (SLOPE^(X-X0)+Y0),COLOR); 
END; (♦FOR^) 
END; (♦ELSE^) 
END (♦IF^) 
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ELSE IF (Y0>Y1) THEN BEGIN 
FOR Y:-YO DOWNTO YI DO BEGIN 

PLOTPOINT (X0,Y,COLOR); 
END (♦FOR*) 
END (♦IF*) 
ELSE BEGIN 
FOR Y>Y1 DOWNTO YO DO BEGIN 

PLOTPOINT (X0.Y.COLOR); 
END (♦FOR*) 
END(*IF THEN ELSE*) 
END; 

BEGIN 

WRrrELN(OUTPUT,'ENTER THE BEGINNING TIME.'); 
WRITELN(OUTPUT,' IT MUST BE A REAL NUMBER GREATER THAN', 
•OR EQUAL TO ZERO.'); 

READLN (T2); 

WRITELN(OUTPUT,'ENTER THE TOTAL TIME THAT THE OBJECT, 
• FALLS.'); 

WRrrELN(OUTPUT,'THIS VALUE MUST BE GREATER THAN THE", 

• BEGINNING TIME.'): 
READLN (LAST); 

WRITELN(OUTPUT,'ENTER THE INTERVAL YOU WANT THE TIME', 

' TO BE DIVIDED BY.'); 
WRITELN(OUTPUT,'USE A DECIMAL TO REPRESENT A FRACTIONAL', 

' PART OF A SECOND.'); 
READLN (DELTAT); 

WRITELN(OUTPUT,'TIME DISTANCE'); 
T:-T2; 

WHILE (T<-LAST) DO BEGIN 
S>(G*T*T)/2; 

WRITELN(OUTPUT,T: 10:4,S: 10:4); 

T:-T+DELTAT; 
END;(*WmLE*) 
READLN; 

WRITELN(OUTPUT,'IF YOU WOULD PLEASE WAIT I ', 

WILL TRY TO DRAW A GRAPH OF THESE VALUES.'): 
DELAY(2000); 
T:-T2; 

GRAPHICSON; 
CLEARGRAPHICS; 
PRCDRW(0,0,XMAX.0, 1); 
PRCDRW(0,0,0,YMAX,1); 
WHILE(T<-LAST) DO BEGIN 
S:-(G*T*T)/2; 

PLOTGENERAL (T,S,0.0,LAST0.0,(G*LAST*LAST)/2, 1); 
T>T+DELTAT; 
END;(*WHILE*) 
READLN; 
GRAPHICSOFF; 
END. 
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Sample Output 
Running 

ENTER THE BEGINNING TIME. 

IT MUST BE A REAL NUMBER GREATER THAN OR EQUAL TO ZERO 
0.0 

ENTER THE TOTAL TIME THAT THE OBJECT FALLS. 

THIS VALUE MUST BE GREATER THAN THE BEGINNING TIME. 

1.0 

ENTER THE INTERVAL YOU WANT THE TIME TO BE DIVIDED BY. 
USE A DECIMAL TO REPRESENT A FRACTIONAL PART OF A SECOND. 
0.01 

TIME DISTANCE 
0.0000 0.0000 
U.0100 0.0005 
0.0200 0.0020 
0.0300 0.0044 
0.0400 0.0078 
0.0500 0.0123 
0.0600 0.0176 
0.0700 0.0240 
0.0800 0.0314 
0.0900 0.C397 
0.1000 0.0490 



0.9000 3.9690 

0.9100 4.0577 

0.9200 4.1474 

0.9300 4.2380 

0.9400 4.3296 

0.9500 4.4223 

0.9600 4.5158 

0.9700 4.6104 

0.9800 4.7060 

0.9900 4.8025 
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NUMERICAL PROGRAMS 

NUMERICAL PROGRAM 1: Monte Carlo Area 
AUTHOR: Randy Odendahl 

OBJECT To use a Monte Carlo method to estimate the area of a circle. This is done 
by using the monitor screen (size 3:^0 by 200 « 64000 square pixels) and a circle of 
radius 100 on the screen, x is chosen uniformly between 0 and 319 and y is chosen 
uniformly between 0 and 199. Then if we repeat this count 20000 times we see that 

Area of circle _ Dots in circle _ Dots in circle 
Area of screen " TotalDots " 20000 

. r • I Dots in circle ^ . - 

Area of circle = ♦ Area of screen 

20000 

A^^^ - I - Dots in circle ^ 

Area of circle = — — ♦ 64000 

20000 

Therefore, Area of circle = 32 ♦ Dots in circle 



Listing of monteCarloArea: 

program monteCarloArea; 
{monteCarloArea uses Monte Carlo methods 
to estimate the area of a circle.} 

{$Ib:graphics.pas} {See appendix for a listing of graphics.pas} 

const 

xcenter* 320; 

>center» 100; 

radius- 100; 

totalCount- 20000; 
var 

x»y: real; 
insideCount: integer; 
i: integer; 

begin 

randomize; clearGraphics; graphicsOn; 
insideCount:* 0; 

for i> 1 to totalCount do begin 
x> random(xmax); 
y> random(ymax); 
if (( sqr(x-xcenter)+sqr(y-ycenter)) 
<» sqr(radius) ) then begin 
insideCount:« insideCount +1; 
plotPoint(round(x)^und(y), 1); 
end {if}; 
end {for}; 
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dday(lOOOO); 
graphicsOfr; 

writeln(output, 'There were MnsideCountl, 

' points inside the circle and'); 
writeln(ouQ)utt' '.totalCoum-insideCount): 1, 

' points outside'); 
writeln(outputt'The estimated area of the circle', 
' is ',( (1.0*insideCount/totalCount)* 
(I.O*xmax*ymax) ):I:2); 
writeln(output,'The actual area is ',pi*sqr(radius):l:2); 
end. 



Listing of sample output: 

There were 10Q28 points inside the circle 
and 9972 points outside. 
The estimated area of the circle is 32089.60 
The actual area is 31415.93 
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NUMERICAL PROGRAM 2: Minimize area. 

SOURCE: Dwayne Channel and Christian Hirsch, "Computr Methods for Problem 
Solving in Secondary School Mathematics". 1984 NCTM Yearbook, pp. 174-175. 
(See Bibliography for a full reference for the Yearbook.) 

Translated into Pascal by Randy Odendahl. 

OBJECT To investigate the opdmal radius and height of a cylindrical container of 
fixed volume to minimize its surface area. 

V ^ KrH 

Write h and Surface as functions of V and r. 

S = 2jcrA +2jcr* = 2;cr(A + r) 

Listing of minimlzeArea: 

program minimizeArea(input,output); 
const pi-i3.141S926; 

var radius, height, surfaceArea, deltaRadius, volume: real; 
totalTrials, trial: integer; 

begin {minimizeArea} 
volume:«236.0; radius:»0.5; deltaRadius:»0.5; totalTrials:»20; 
writeln(output, 'Program to investigate c^timal dimensions to minimize '); 
writeln(ou^t,'surface area of a cylindrical container of fixed volume '); 
writeln; 

writeln(output,'Enter fixe^' volume of the cylinder.'); 
writeln(ouq)ut,volume:5: 1); 
writeln(output,'Enter initial length of radius.'); 
writeln(ouq)ut»radius:3: 1); 
writeb(output,'Enter increment for radius.'); 
writeln(ouq)ut,deltaRadius:3: 1); 
writeln(output,'Enter number of trials to investigate.'); 
writeln(ouq)ut,totalTrials); 

writeln(output,'VOLUME - '.volume: 1:2/ TRIALS .'.totalTrials:!); 
writeln(outpui,' RADIUS HEIGHT SURFACE AREA'); 
writeln(o»iq)ut,' - '); 

for trial:* 1 to totalTrials do begin 
height:* volume/(pi*sqr(radius)); 
surfaceArea:* 2*pi*radius*(height + radius); 
writeln(ou!put.radius:10:2,height:10:2,surfaceArea:lS:2); 
radius:* radius +deltaRadius; 
end {for} 
end. {minimizeArea} 
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Sample output; 

Program to investigate optimal dimensions to minimize 
surface area of a cylindrical container of fixed volume. 

Enter fixed volume of the cylinder. 
236.0 

Enter initial length of radius. 
0.5 

Enter increment for radius. 
1.0 

Enisr number of trials to investigate. 
10 

VOLUME - 236.00 TRIALS -20 
RADIUS HEIGHT SURFACE AREA 



0.50 


300.48 


945.57 


1.50 


33.39 


328.80 


2.50 


1102 


228.07 


3.50 


6.13 


211.83 


4.50 


3.71 


232.12 


5.50 


2.48 


275.88 


6.50 


1.78 


338.08 


7.50 


1.34 


416.36 


8.50 


1.04 


509.49 


9.50 


0.83 


616.74 



NUMERICAL PROGRAM 3: Bisection 
AUTHOR: Lynn R. Ziegler 

OBJIflCT The following program is a graphics picgram to graphically show how the 
method of bisection can be used to find a zero of a function. In our case the function 
is y=2-e' in die range Q^^l. We expect to find the root at ln(2)-0.6931472... . 

Bisection works by evaluating die function at die ends of ai interval where die func- 
tion has different signs (positive at one end, negative at the odier). The function is 
dien evaluated at the mi(4x>int of that interval. If die function is zero die root is found; 
if positive, die new interval will be between the midpoint anJ d^^ end having a nega- 
tive functional value; if negative, the new interval will be between die midpoint and 
die end having a positive functional value. This is continued, shrinking the interval of 
interest until it becomes smaller dian some tolerance. Then die root is at the middle 
of die small interval remaining plus or minus half diat tolerance. 

This program does ^he bisection graphically, showing die intervals cm screen as higher 
and higher "walls'* closing in on die root 
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Listing of Program Bisection 
program bisection(input,output); 

{$igraphics.pas} {These two statement include some grs^hic^ 
{$Igra(^ics.cwo} {needed to nm the program in Turbo Pascal on an IBM PC) 
{Listings of them can be found in the appendix.} 

const xmin-0.0; 
ymin—2.0; 
xlarge-LO; 
yliirge-2.0; 

var x4ower,upper,mid,dot:real; 
signk)w,i: integer; 

function f(x:real):real; 

{This function will be used in the bisection to 
compute valuta for finding the solution of f(x)-iO} 

begin {f} 

f:-2.0-exp(x) 
end; {f} 

function sign(y:real):integer; 
begin {sign} 

if (y-0.0) then sign>0 

else if (y<0.0) then sign:— I 

else sign>I 
end; {sign} 

begin {bisection} 
graphicson; 

plotaxe$(xmin»xlai;ge»ymin»ylarge, I): 
x:«xmin; 

while(x<«xlarge) do begin 
plotgeneral(x»f(x)»xmin,xlarge, 

ymin.ylarge,l); 

x>x-h0.0016 
end; 

dot>0.0; lowen-0.0; upper> LO; 
signlow>sign(f(0.0)); 
while((upper-lower) > 0.01) do be^in 
mid: «(upper+lower)/2.0; 
dot:-dot-i0.15; 
if sign(f(mid))«0 then begin 
lower-mid; 
uppen*mid; 

drawlinegeneral(xmin,dot,lower»dot,xmin, 

xlarge,ymin,ylarge,l); 
drawlin^eneral(upper,0.0»upper,dot,xmin» 

xlarge,ymin,ylarge.l) 

end {iO 
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else if (sign(f(inid))»signlow) dien begin 
lowen»mid; 

diawlinegeneral(xnun,dot,lower,doumin» 

xlarge,yniin,ylarge,I); 

drawlinegeneral(lower,O.OJower,doumin, 

xlarge,ymin,ylarge, 1) 

end 
else begin 
uppen-mid; 

diawlinegeneral(upper,dot,xla]:ge,dot,xmin, 

xlarge,ymin,ylarge,I); 

drawlinegeneral(upper,0.0,upper,doumin, 

xlaige,ymin,ylarge,l) 

end 

end; {w!iue} 

writelnCThe lower x,y pair is: \ 

lowcr8:3,f(lower):8:3); 
writelnCThe upper x,y paii is: \ 

upper.8:3,f(upper):8:3); 

delay(SOOO); 
graphicso/F 
end {bisect} 



Sample Output 
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Appendix 



Graphics Procedures 
AUTHOR: LynnZieglcr 



This sq)pendix contains standard types and procedures u.«cd by graphics programs written for use with 
Turbo Pascal on the IBM PC. Similar prognuns were written for use with Turbo Pascal on Apple II 
computers with CP/M boards. Copies of either the IBM or Apple routines are available on request 



The following are die constants, types, and routines in the file "GRAPHICS.PAS" which is included in 
many of die above samjde programs. 

const 

xmax-319; 
ymax-199; 
pi - 3.1415926535; 



type 

xvalue-0..xmax; 
yvalue-0..ymax; 

functionarray«array[0..xmax] of yvalue; 



proceoure cleargraphics; 
begin {cleargraphics} 

graphcolormode; 
end; {cleargraphics} 



procedure graphicson; 

{This procedure prepares tht system for use of graphics.} 

begin {graphicson} 
gn^hcolormode; 
end; {graphicson} 



procedure gnq)hicsoff; 

{This procedure restores normal text mode.} 

begin {graphicsofT} 
Textmode(bw80) 
end; {graphicsoff} 
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|m)cediire plotpoint(x,y,cdor : integer); 

{This procedure plots a point at coordinates x,y on the screen. If color 
is 0 the tiackground will be plotted. Otherwise white will be plotted.} 

var xtenip t xvalue; 
ytemp : yvalue; 

begin {ploipoint} 

if (((k-x) and (x<-xmax) and (0<-y) and (y<-ymax)) then begin 
xtenip>abs(x) mod (xmax+1); 
ytenip:-ymax - (abs(y) mod (ymax+1)); 
pk)t(xtemp»ytemp,color) 
end {if} 
end; {plotpoint} 



procedure smoothplot(f : funcdonarray; coIcmt : integer); 
{This procedure accepts an array f of xmax int^ers in the range 0 to ymax and 
plots them on the screen. It assumes the integers rei»esent a continuous 
function so it smooths vertical jumps as much as it can.} 

var ij»mid : integer; 

begin {smoothploc} 
ploi(04I0],color); 
for i:«l to xmax do begin 

{Is there a junq> where the current value is more than one dot 
above the previous value? If so, plot the vertical values to 
fill in between the two points.} 

if f[t]>f[i-l]+l then begin 
mid:-(f[t-l]+f[i]) div 2; 
fr/ j:«f[t-l] to mid do plotpointvi-lJ,color); 
tor j>f[t] downto (mid+l) do ploq)oint(ij,color) 
end {if} 

{Perh^ there is a jump down at least two dots? If so, fill in 
the vertical values as above.} 
else if f[i]<f[i-l]-l then begin 
mid:-(fli-l]+fli]) div 2; 

for j:«f[i-l] downto (mid+1) do plo^int(i-lJ,color); 
for j:»f[i] to mid do plotpoint(iJ,color) 
end {elseif} 

{Maybe no jump occun. Then just plot the point} 
else plotpoint(i,f[i],color) 

end {for} 
end; {smoothplot} 
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The following are the constants, types, and routines in the file ''GRAPHICS.TWO'* which is included in 
many of the above sample programs. 

type 

arrayofTeals-array[1..120] of real; 



imcedure drawUne(xO,yO»xl,yl, color : integer); 

{This procedure draws a line between the points (xO.yO) and (xl,yl) with shade 
color. The screen has coordinates 0 to ymax from bottom of screen to top of 
screen and 0 to xmax from left to right The drawUne uses full screen 
density and plots nice dense lines by checking slopes.} 

var 

x,y : integer; 
slope : real; 

begin {drawline} 

if (xO-xl) then {Vertical lines are plotted by simply changing y 
and plotting the same x value (xO) at all points.} 
if (yO<yl) then 

fof y>yO to yl do plotpoint(xO,y .color) 
else 

for y>yO downto yl do plotpoint(xO,y,color) 
else begin {Non vertical lines are plotted here.} 
slope:-(yl-yO)/(xl-xO); 

if (abs(slope)<-1.0) then {Small slopes are handled by varying 
X and computing appropriate y values.} 

if (xCkxl) then 
for x>xO to xl do 

plotpoint(x»round(slope*(x-xO)+yO)»color) 

else 

for x>xO downto xl do 

plotpoint(xjound(slope*(x-xO)+yO)»color) 
else begin {Laige slopes are handled by thinking of x 

as a function of y, varying y» and then 
computing the proper x values.} 

slopei-l.O/slope; 
if (yO<yl) then 
for y:-yO to yl do 

plotpoint(round(slope*(y-yO)+xO),y,color) 

else 

for y>yO downto yl do 

plotpoint(round(slq)e*(y-yO)+xO),y,color) 
end {if-then-else} 
end {if-then-else} 
end; {diawline} 
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procedure plotgeneral(x»y»xinin.xlaige»yinin»ylai8e : real; color : integer); 
{This procedure plots a point at position (x»y) in the coordinate system which 
places xnrJn as the leftmost x value* xlarge as the rightmost x value* ymin as 
the bottom y value* and ylarge as die top y value. It is plotted in shade 
color. It uses proportiops and die specified values for standard screen 
coordinates (te., stand'ird coordinates are x values 0 to xmax from left to 
right and 0 to ymax fr:m bottom to top).} 

begin {{dotgeneral} 

ploq)oint(n>und((x-xmin)/(xlarge-xmin)*xmax), 

round((y-ymin)/(ylaige-ymin)*ymax), color) 
end; {plotgeneral} 



procedure drawlinegeneraI(xO,yO,xl,yl,xmin^large,ymin»ylarge : real; 

color : integer); 

{drawlinegeneral (Laws a line between (xO,yO) and (xl»yl) in the coordinate 
system running from xmin at left to xlarge at right and ymin at bottom and 
ylarge at top. It works by converting the points (xO,yO) and (xl»yl) to 
standard screen coordinates and then calling jmcedure drawline.} 
var 

xfactCM*»yfactor : real; 

begin {drawlinegeneral} 

xfactCMr:»xmax/(xlarge-xmin); 
yfactKM:»ymax/(ylarge-ymin); 

drawline(round((xO-xmin)*jtfactor)»round((y()-ymin)*yf?jto^^ 

round((xl-xmin)*xfactor)4ound((yl-ymin)*yfactor)»color) 
end; {drawlinegeneral} 



procedure drawcirclegeneral(xcenter»ycenter^adius : real; 

xmin»xlarge,ymm»ylarge : real; color : integer); 

{This procedure draws a circle centered at (xcenter»ycenter) of radius radius. 
Tlie circle is in relative coordinate system xmin left to xlarge right and ymin 
bottom to ylarge top. The circle is drawn by using polar coordinates centered 
at (xcenter»ycenter). The relative number of points is proportional to the 
'^normalized'* radius and the number of points on dis screens perimeter.} 

var 

theta,deltatheta : real; 

begin {drawcirclegeneral} 
theti^-O.O; 

deltatheta>2.0*(xlarge-xmin+ylarge-ymin)/radius/(xmax+ymax); 
while (theta<»2.0*pi) do begin 

plotgeneral(xcenter+radius*cos(theta)»ycenter+radius*sin(theta)» 

xmin»xlarge»ymin»ylarge»color); 
theu:»theta+deltatheta 
end {while} 
end; {drawcirclegeneral} 
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procedure plotaxes(xmin,xlarge,ymin,ylarge:real; colorinteger); 

{Tills procedure draws x and y axes through the origin (0,0) in the coordinate 
system defined by xmin left to xlarge right» ymin bottom to ylarge top in 
shade color. It converts to standard screen coordinates to draw the line.} 

begin {plotaxes} 

drawIin^eneral(0.0,ymin,0.0,ylaige,xnun,xlarge,ymin,ylarge,color); 
drawlinegeneral(xmin,0.0,xlarge,0.0»xniin»xlarge,ynun,ylarge,color) 
end; {plotaxes} 



procedure ConvertToPblar(x,y.xcenter,ycenterreal; var radius,theta:ieal); 

{This procedure converts the point (x,y) to its represenutioa in polar 
coordinates via the pair (radius^theu). The polar coordinates are with 
reference to axes through (xcenter,ycenter). Standard transformations are 
used for conq)uting both radius and theta.} 

begin {ConvertToPdlar} 

radius>sqrt(sqr(x-xcenter)+sqr(y-ycenier)); 
if x-xcenter then 
if y>ycenter then theta:-pi/2 
else theta:-l.S*pi 
else if x>xcenter then theta:-arctan((y-ycenter)/(x-xcenter)) 
else theta:-arctan((y-ycenter)/(x-xcenter))+pi 
end; {ConvertToPolar} 



procedure Rotate(var x,y:arrayofTeals; theta:real; 

numberofjpointsiinteger; xcenter.ycenterreal); 

{This procedure causes the numberofjpoints points (x[i],y[i]) in the arrays x 
and y to be rotated by angle theta with respect to axes intersecting at 
(xcenter,ycenter). Contents of arrays x and y wUl be changed.} 

var 

i : integer; 
radius,thetaO : real; 

begin {Rotate} 

for i:-l to nuniberofpoints do begin 

ConvertToPolar(x[i],y[i],xcenter,ycenterjadius,thetaO); 
x[i]:-xcenter+radius*cos(thetaO+theta); 
y[i]>ycjnter+radius*sin(thetaO+iheu) 
end {for} 
end; {Rotate} 
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